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Abstract 

Scintillation efficiency of low-energy nuclear recoils in noble liquids plays a crucial 
role in interpreting results from some direct searches for Weakly Interacting Massive 
Particle (WIMP) dark matter. However, the cause of a reduced scintillation effi- 
ciency relative to electronic recoils in noble liquids remains unclear at the moment. 
We attribute such a reduction of scintillation efficiency to two major mechanisms: 
1) energy loss and 2) scintillation quenching. The former is commonly described 
by Lindhard's theory and the latter by Birk's saturation law. We propose to com- 
bine these two to explain the observed reduction of scintillation yield for nuclear 
recoils in noble liquids. Birk's constants kB for argon, neon and xenon determined 
from experimental data are used to predict noble liquid scintillator's response to 
low-energy nuclear recoils and low-energy electrons. We find that energy loss due 
to nuclear stopping power that contributes little to ionization and excitation is the 
dominant reduction mechanism in scintillation efficiency for nuclear recoils, but that 
significant additional quenching results from the nonlinear response of scintillation 
to the ionization density. 
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1 Introduction 



Noble liquid scintillators such as liquid xenon pQ , argon [2|3] , and neon [I] are 
expected to be excellent targets and detectors for direct dark matter detec- 
tion experiments searching for Weakly Interacting Massive Particles (WIMPs), 
which may constitute the dark matter in the universe [5TI61T71I8] . These exper- 
iments measure scintillation light induced by low-energy nuclear recoils due 
to elastic scattering of WIMPs. Absorption of nuclear recoil energy in noble 
liquid scintillators produces excitons and electron-ion pairs along the track. 
Free excitons collide with ground states to form excited molecules (excimers) 
through 

X* + X -»■ X*, (1) 



where X stands for any type of noble liquid. Free ions undergo collision, re- 
combination and deexcitation processes, 

Xt + e~ -> X** + X, 

(2) 

X** -> X* + heat, 



X* + X -> X, 



2 ■ 



to form excimers. The excimers then decay radiatively from the lowest-excited 
molecular states 1 £+ and 3 £+ to the repulsive ground state X E+. 

It is well known that noble liquid scintillators have reduced scintillation yield 
for low-energy nuclear recoils compared to electronic recoils [9)10)11)12)13)14] . 
Only a fraction of the energy loss results in ionization and atomic excitation. 
Moreover, high ionization density undermines recombination of electron-ion 
pairs and reduces scintillation light yield. The relative scintillation yield, de- 
fined as the ratio of the numbers of photons emitted from pure noble liquids in 
nuclear and electronic recoil events at the same energy, is a good measurement 
of the nuclear recoil scintillation efficiency, g/, determined by the visible de- 
posited energy over the true recoil energy. In the case of electrons and 7-rays, 
almost all the energy loss by ionization is converted into scintillation light 
through electron-ion recombination, so the relative scintillation efficiency is 
assumed approximately equal to 1 in the absence of an electric field. But in 
the case of nuclear recoils, qj is much smaller than 1 and can vary as a function 
of nuclear recoil energy. 

The scintillation efficiency plays an important role in the direct detection of 
WIMPs. In the design of a new experiment, the scintillation efficiency is related 
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to the detection threshold and hence to the background level and ultimate 
sensitivity. In the interpretation of an experimental result, the scintillation 
efficiency is crucial to the determination of WIMP mass and WIMP-nucleon 
cross section. The nuclear recoil scintillation efficiencies for liquid xenon, ar- 
gon, and neon have been measured [9lll0|lll|ll2|ll3|ll4j using neutron sources. 
The detector is usually calibrated with well-known 7-rays, such as 122 keV 
and 133 keV lines from a 57 Co source. The relative scintillation efficiency for 
7-rays, e 7 , is defined as the visible energy divided by the incident 7-ray en- 
ergy, and is assumed to be 1. As discussed later in this paper, this assumption 
is valid for electronic recoil energy above 20 keV in the noble liquids under 
investigation. The nuclear recoil scintillation efficiency is determined by the 
ratio of the nuclear recoil visible energy using the electron-equivalent energy 
calibration to the true recoil energy E R , qf = E R m /E R . 

The measurements are usually compared to either Lindhard's theory [T3] or 
Hitachi's treatment [16]. It was found that Lindhard's theory alone can not 
well explain the observed behavior in the data. The alternative explanation, 
Hitachi's treatment, states that a biexcitonic quenching mechanism can oc- 
cur before the free excitons self-trap when the excitation density is very high. 
This explanation can agree reasonably well with the xenon data [9], but has 
not been applied to explain low-energy recoil data for neon or argon. More- 
over, there are other possible quenching processes that could contribute to 
the reduction of scintillation efficiency for low-energy nuclear recoils. These 
include collisions (via the Penning process |17j ) between two excited molecu- 
lar states (excimers) to form one excited state and one ground state [IB] , and 
supereleastic collisions that quench the singlet states to triplet states |19j . 

A more universal description of the reduced scintillation efficiency for nuclear 
recoils is preferred for all noble liquid scintillators. Since the proposed quench- 
ing mechanisms are all dependent on the density of the ionization and excita- 
tion track left by the recoiling nucleus, it is possible to form a combined model 
of these mechanisms without incorporating details of the relative contributions 
of the different mechanisms. Birk's saturation law for organic scintillators [20] 
provides a convenient description of the dependence of scintillation quench- 
ing on ionization density. In this study, we apply Birk's law to noble liquids, 
offering a conventional way to determine the total scintillation efficiency of 
nuclear recoils by measuring Birk's constant (kB). It has been shown that 
the luminescence intensity in the noble gas scintillator depends solely on the 
energy density and is independent of the kind of the particle [2~Tf22] . This is 
to say that a measurement of kB for noble liquids will allow understanding 
of the relative scintillation efficiency for nuclear recoils induced by neutrons, 
alphas, and other heavy isotopes. This is a very valuable way to determine 
the scintillation efficiency for nuclear recoils induced by all types of particles 
in noble liquid scintillators. Furthermore, this method allows determination 
of the relative scintillation efficiency for nuclear recoils in noble liquids by 
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measuring the constant kB with 7-rays. This is a much easier measurement 
compared to the nuclear recoil measurements with neutron and alpha sources. 

In this paper, we propose a model to combine Lindhard's theory and Birk's 
saturation law to describe the reduction in scintillation efficiency observed 
in noble liquid scintillators. We describe these two reduction mechanisms in 
Section [2] and [3], respectively. The model combining these two reduction mech- 
anisms is presented in Section H] and its predictions are compared to experi- 
mental data in Section [51 The scintillation efficiency (e 7 as a function of recoil 
energy) for very low-energy electrons and 7-rays is discussed briefly in Sec- 
tion O Finally, we summarize our conclusions in Section [7J 



2 Reduced Ionization Energy by Nuclear Collisions 

When a neutron or WIMP scatters elastically off a noble liquid atomic nu- 
cleus, the recoiling nucleus then loses its energy by colliding with electrons and 
nuclei within the detector. This nuclear recoil process involves the competi- 
tion between, on the one hand, energy transfer to atomic electrons and, on the 
other hand, energy transfer to translational motion of atoms. The total rate 
at which the recoiling nucleus loses energy with respect to distance (dE/dx) 
is dependent on the medium through which it travels, and is also called the 
stopping power. At low energies, the total stopping power of the noble liquid 
atom consists of electronic and nuclear stopping power. The electronic stop- 
ping power is the amount of energy per unit distance that the recoiling nucleus 
loses due to electronic excitation and ionization of the surrounding noble liq- 
uid atoms. The nuclear stopping power is the energy loss per unit length due 
to atomic collisions that contribute to the kinetic energy (thermal motion) of 
the noble liquid atoms, but that do not result in internal excitation of atoms. 
The proportion of electronic to nuclear stopping power depends on the recoil 
energy of the nucleus. If the recoil energy were very large, the nuclear stop- 
ping power would be very small compared to the electronic stopping power. 
However, in the energy range of WIMP-nucleus elastic scatterings, the nuclear 
stopping power plays a significant role in the energy loss of the recoiling noble 
liquid nucleus. J. Lindhard et al. [15] discussed in detail the theory of energy 
loss of low-energy nuclei. 

Supposing that the recoiling nucleus loses all of its energy in the detector, the 
total energy loss can be expressed in terms of the losses due to the electronic 
stopping power rj and nuclear stopping power v as [TB] 

E R = r){E R )+v(E R ), (3) 
where rj and v are both functions of recoil energy E R . As only the portion of 
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the energy lost in electronic excitation or ionization will result in the creation 
of excitons and electron-ion pairs in the noble liquids, the fraction defines an 
ionization energy reduction factor (f n ) due to losses to the nuclear stopping 
power 

UEr) = ^n~~ v(En) + u(E R y (4) 



As the total stopping power is 
dE _ dE dE 

I - ! Jtot — {—T~ Jelec + (.—r~ Jnucb (O) 

dx dx dx 



f n {E R ) can then be determined by the ratio of two integrals 

f(E)= J ER (dE/dx) clcc dE 

I ER ((dE/dx) clec +(dE/dx) nud )dE' 



To present f n as a function of recoil energy, the integrals above should be 
evaluated for each possible recoil energy. Lindhard et al. (T5] represents f n as 

kg(e) 

U l + kg(ey [0 



where, for a nucleus of atomic number Z, e = 11.5Er (keV) Z 7 / 3 , k = 
O.mZ 2 / 3 ^ 1 / 2 , and g(s) is well fitted by: g(e) = 3e a15 + 0.7e - 6 + s. Fig. Q] 
shows this ionization energy reduction factor for noble liquids from Lindhard's 
theory. 



3 Reduced Scintillation Yield due to High Ionization Density 



3.1 Birk's Saturation Law 



The passage of a particle in a noble liquid produces a structured track along its 
path that is conveniently described in terms of a core and a penumbra [23] . The 
penumbra that surrounds the core is a low ionization density zone. The core is 
expected to be a high ionization density zone, so ionization density dependent 
quenching caused by biexcitonic collisions or the Penning process is likely to 
occur there. On the one hand, the free excitons can be self-trapped to form 
excimers, for example, Ar%, that then fluoresce to the lowest states. On the 
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Fig. 1. The ionization energy reduction factor according to Lindhard's theory versus 
the recoil energy in liquid neon (black line), argon (red dotted line) and xenon 
(magenta dashed line). 

other hand, the free excitons can diffuse and undergo biexcitonic quenching, or 
the excimers can further collide with each other via the Penning process, with 
a probability of occurrence that depends on the density of free excitons pro- 
duced, which is proportional to ionization density. Therefore, the biexcitonic 
collisions and the Penning process act as quenching agents for the excitons pro- 
duced by the ionization along the track. It is well known that the number of 
excitons and electron-ion pairs produced per unit path length is proportional 
to the electronic energy loss dE/dx [23], with a proportionality constant that 
will be designated as A. The local concentration of the core is also propor- 
tional to the ionization density, and is given by BdE/dx. The overall collision 
probability in the core is given by k. Thus, the specific fluorescence can be 
expressed by 

dS = Af 

dx 1 + kB^' 1 ; 

dx 



Eq. [8] is called Birk's saturation law [20] , which describes the relative scintilla- 
tion response of scintillators to an ionizing particle of any energy. The values 
of A and kB can be determined experimentally. 

According to Eq. [HJ the scintillation photon yield is reduced at high ionization 
density. Thus, we define a quenching factor 



fl i , kB dEi (9) 



dx 



6 



which is related to the electronic stopping power dE/dx induced by nuclear 
recoils. 



3.2 Electronic Stopping Power for Heavy Ions 



For ions heavier than protons, the electronic stopping power in a given material 
can be calculated from the stopping power in the same material based on the 
"heavy ion scaling rule" , 

s A = (c^oX (io) 



where the stopping power S p for protons is determined at the same velocity 
as the stopping power Sa for heavy ions. (Z\ is the 'effective charge' for ions 
of atomic number Z\. Effective ion charges will result from the stripping of 
bound electrons from the ion when moving through a medium. Thus, the 
effective charge fraction ( is expected to be related to the relative velocity, 
IV, between the ion velocity V\ and the velocity of the valence electrons in the 
medium. When the valence electron gas of the medium is characterized by 
the effective number of electrons that participate in plasma excitations, the 
velocity of the valence electrons can be characterized by the Fermi velocity 
vf = (37r 2 n) 1 / 3 /i/m, where n is the electron density and m the electron mass. 
v r is proposed in [25] to be 

v r = < ^ 15 ^ h (11) 

[^(l + i^), iivi>v F . 



A formula for the effective charge fraction ( has been deduced by Brand and 
Kitagawa [26] as a good approximation in the region where stopping power is 
proportional to the velocity v\\ 

C . ? + 0.5(1-,) In + (12, 



where vq and a>o are the Bohr velocity and radius, respectively, q denotes the 
degree of ionization, which was calculated by applying a velocity-stripping 
criterion to the Thomas- Fermi model [27] of the neutral atom and is tabulated 
as a function of the reduced variable y r = v T j{Z\ Vq) in [26]. An alternative 
energy stripping criterion was proposed by Mathar and Posselt (MP) [28] 
to explain the dependence of the ionization fraction on the ion velocity and 
to yield better results especially for higher ion velocities. Using the energy 
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definition of the Brandt-Kitagawa (BK) model, q is determined by solving 
numerically the following equation 



80(1 - 1?'V = *2±^, (13) 



where a = 0.24005 and y = v 1 /(v G Z 2 ^). The screening radius A, which is used 
to parameterize the charge density profile of the projectile ion, is determined 
by minimizing the total energy of the bound electrons in BK to be 

A = 2fl ( 1 -^ /3a ° (14) 

^i 1/3 [i-(i-g)/7] 



Ziegler, Biersack and Littmark (ZBL) [29j modified the ion size parameter A 
to include a tabulated factor individual to Z\. The ionization fraction q is then 
parameterized as a universal fitting function of y r as 

1 (0.803y°- 3 -1.3167^ 6 -0.38157y r -0.008983y2) /, -s 



The ionization fractions from the three different approaches described above 
are reproduced in Fig. [2] as a function of the variables y or y r . In this paper, 
we adopt the BK formalism to calculate the effective charge of ions and the 
electronic stopping power values, but the alternative approaches are used to 
estimate the theoretical uncertainties. 
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Fig. 2. Ionization fraction as a function of variable y or y r in the BK (black solid 
line), ZBL (blue dotted line) and MP (red dashed line) models. 

We can utilize Eq{10] and the effective ion charge to calculate the electronic 
stopping power, provided that the stopping power values for protons in the 
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same material are known. A table of stopping powers and ranges for pro- 
tons and alpha particles has been produced by a committee of the Interna- 
tional Commission on Radiological Units and Measurements (ICRU) [3U]. In 
the present study, we retrieve the proton stopping powers from the PSTAR 
database [31]. Fig. [3] shows the electronic stopping power for proton projec- 
tiles in Ar target as a function of the proton velocity (3 = v/c. According to 
Lindhard and Scharff [32], the stopping power at low energy is proportional to 
the projectile velocity. This is illustrated by a proportional function (red thick 
line) fit in Fig. [3J We will use this fitting function to estimate the stopping 
power at lower energy, where it is not covered by the databases, to produce 
the electronic stopping power for ions at kinetic energy as low as 1 keV. 
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Fig. 3. Electronic stopping power for protons in an Ar target versus the velocity 
/3 = v/c of the proton projectiles. A proportional function (red thick line) is fitted 
to the low velocity region with a parameter of 5.776 x 10 4 ± 349.7 MeV cm 2 /g. 

Fig. H] gives the electronic stopping power for Ar ions in Ar targets. The red 
solid line is produced based on BK model and the "heavy ion scaling rule". 
Paul and Schinner (PS) [33] performed a systematic comparison of a large 
collection of stopping power data for projectiles from 3 Li to 18 Ar to those for 
alpha particles in the same materials. They found that, for each target element, 
a slightly adapted sigmoid function with three parameters can describe the 
normalized relative stopping power 



rcl 



s A /zl 

Shc/2 2 ' 



(16) 



For comparison, results from the PS empirical scaling approach have also 
been plotted in Fig. [4] with a blue dashed line. The difference between the two 
calculations is about 10% at kinetic energy of 1 MeV and increases to 35% at 
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10 keV. In addition, a green dash-dotted line in the figure shows a calculation 
based on the "SRIM" program of Ziegler [29]. The results deviate from those 
of the BK and PS approaches at higher energy with a kink at ~ 0.1 MeV, but 
are comparable to the results of the BK model at lower energy. 



Fig. 4. Electronic stopping powers for Ar ions in Ar targets as a function of the 
kinetic energy. The red solid line is from the BK model, the blue dashed line from 
the PS empirical scaling approach and the green dash-dotted line from Ziegler's 
SRIM program. 

Based on the BK model, we show in Fig. [5] the electronic stopping power for 
Ar (black solid line), Ne (blue dotted line) and Xe (red dashed line) ions with 
themselves as the target materials. By comparing different effective ion charge 
calculation scenarios, the theoretical uncertainties on the stopping power val- 
ues are estimated to be about 15% at 1 MeV and increase to about 50% at 1 
keV where the theoretical description is known to become less reliable. 



4 Model Prediction 

The fraction of nuclear recoil energy, f n (see Eq. [7j), described in Section 
[2] ionizes or excites atoms in noble liquids. Ionizing particles produce exci- 
tons, electron-ion pairs, and a localized concentration core along the track. 
There are two cases: 1) free excitons form excimers through a self-trapping 
process [3"4l35|l36j . and 2) the ions are localized through the formation of ex- 
cited molecular ions and eventually form excimers through recombination, 
deexcitation, and collision. The origins of the luminescence for both cases are 
attributed [371138] to low excited molecular states, namely, or 3 S+. The 
fraction of excitons that undergo self-trapping depends largely on ionization 
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Fig. 5. The electronic stopping powers based on the BK model for Ar (black solid 
line), Ne (blue dotted line) and Xe (red dashed line) ions with themselves as the 
target materials. 

density. The fraction of ions that combine with ground states to form excited 
molecular ions is proportional to ionization density. Furthermore, the fraction 
of excimers relaxing to the lowest-excited states is also a function of ioniza- 
tion density. This is to say that some of the excitons or ions can collide each 
other without going to excimers and some of the formed excimers can further 
undergo the Penning process to reduce the number of radiative lowest-excited 
states. This part of excitons, ions, and the formed excimers is quenched. The 
quenching factor, f h is expressed in Eq. [9j Since f n and f\ are independent of 
each other, the total scintillation efficiency in noble liquids can be represented 



We calculated f n , the fraction of nuclear recoil energy that contributes to the 
ionization process in noble liquids, as a function of nuclear recoil energy, as 
shown in Fig.[U The fraction of the collided excitons or ions, kB, is determined 
from experimental data based on Birk's saturation law by using Eq. [9j 

• For liquid argon, kB = 7.4 x 10 -4 MeV -1 g cm -2 . This is determined by 
using a quenching factor (46%) from a heavy ion measurement [39J, and 
dE/dx (1586.4 MeV cm 2 /g corresponds to 31.9 MeV/amu) is calculated by 
this work. 

• For liquid neon, the Birk's constant, kB = 1.12 x 10~ 3 MeV" 1 g cm -2 , 
is determined by using a quenching factor measured in Ref. [H] with the 
dE/dx values calculated by this work. 



by 



Qf = fn X fl- 



(17) 
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• For liquid xenon, the quenching factor has been measured by several groups [91TT0 
We use the measured quenching factor corresponding to a 70-keV recoil en- 
ergy in Ref. [10] to determine the Birk's constant kB to be 2.015 xlO -3 
MeV -1 g cm" 2 . 

The quenching factor due to scintillation quenching, fi, is then calculated 
according to Eq. [9] for argon, neon and xenon. Multiplying these two reduction 
factors results in the total scintillation efficiency (Eq. [TTj) as shown in Fig. [61 




200 250 300 

Recoil Energy (keV) 



Fig. 6. The total scintillation efficiency for nuclear recoils as a function of recoil 
energy in liquid neon (black line), argon (red dotted line), and xenon (magenta 
dashed line). 



5 Comparison of model prediction with data 



The total scintillation efficiency for each noble liquid, argon, neon, and xenon, 
is compared to experimental data in Fig. [71 Fig. [HJ and Fig. [9j As can be 
seen, the experimental data can not be explained by Lindhard's theory alone, 
but can be well described when the Birk's saturation effect has been taken 
into account. Note that the scintillation efficiencies for 5.035 MeV a particles 
and 33.5 MeV/amu 18 ions in liquid argon were measured to be 71% and 
59% [10139] . which agree very well with our model prediction, 72% and 63%. 

It is worthwhile to mention that, although the theoretical uncertainties of elec- 
tronic stopping power values at very low recoil energy are considerable, the 
effect on the uncertainties of the final quenching factor are not large as the 
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Fig. 7. Argon: a comparison of the total scintillation efficiency (blue curve) cal- 
culated by combining two reduction factors to the experimental data [T2] in liquid 
argon. Also shown as the black curve is the Lindhard ionization energy reduction 
factor. 
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Fig. 8. Neon: the total scintillation efficiency (blue curve) calculated by taking into 
account two reduction factors is presented in liquid neon. 

Birk's constants kB are quite small. For example, by assuming a factor of two 
difference in the electronic stopping powers at 1 keV in liquid argon, the con- 
tribution to the uncertainty of the quenching factor is less than 7%. Thus, the 
dominant uncertainties come from the uncertainties in the Lindhard theory 
at very low energy. However, it is usually believed that the Lindhard the- 
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Fig. 9. Xenon: a comparison of the total scintillation efficiency (blue curve) cal- 
culated by combining two reduction factors to the experimental data |9|10|llj in 
liquid xenon. Also shown as the black curve is the Lindhard reduction factor. 



ory describes experimental data well [4"T|42f 43ll4"4"] , although the experimental 
uncertainties are quite large at low energy. 



6 Low-Energy Response to Electrons and 7-rays 



Noble liquids' response to low-energy electrons and 7-rays is usually assumed 
to be linear for the existing experiments [Tf5f4] . To investigate whether such 
an assumption holds, we also calculate and show in Fig. [10] the quenching 
factor for electronic recoils in liquid argon, neon and xenon utilizing the kB 
constants determined above. It is clear that the assumption holds quite well 
down to 20 keV because of the small stopping power values (a few MeV cm 2 /g) 
above 20 keV. However, the quenching effect is expected to be significant when 
the electronic recoil energy is below 20 keV due to the larger stopping power 
(more than 10 MeV cm 2 /g). 

Note that the kB constants are determined by Birk's law in which the relative 
scintillation efficiency of nuclear recoils is used, under the assumption that the 
electronic recoil quenching factor is 1. As we can see in Fig. [TQl this assumption 
is only valid above 20 keV. 
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Recoil Energy (keV) 



Fig. 10. The quenching factor calculated for electronic recoils in liquid argon (red 
dotted line), neon (black line), and xenon (magenta dashed line). 



7 Conclusion 



In summary, we found that the reduced scintillation response to low-energy 
nuclear recoils comes from two different mechanisms: 1) reduced ionization 
and excitation energy by nuclear collisions and 2) reduced scintillation photon 
yield due to high ionization and excitation density induced by nuclear recoils. 
The former is well described by Lindhard's theory and the latter is attributed 
to biexcitonic collisions between excitons and the Penning process between 
excimers. The scintillation quenching induced by the combination of biexci- 
tonic collisions and the Penning process can be described by Birk's law. We 
have combined these two reduction mechanisms to describe the total scintil- 
lation efficiency for noble liquids. The calculations are compared to available 
data and it is found that they are in good agreement within experimental 
uncertainties. We have also studied the scintillation response of noble liquid 
scintillators to low-energy electrons and 7-rays and found that the response is 
not linear at recoil energy below 20 keV. This paper provides a conventional 
way to measure the total scintillation efficiency of noble liquids for all types of 
particles including neutron-induced nuclear recoils, alpha particles and other 
heavier elements by measuring Birk's constant kB using 7-ray sources. 
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